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ABSTRACT 

We present measurements of positions and relative proper motions in the 30 Doradus region of the 
Large Magellanic Cloud (LMC). We detail the construction of a single-epoch astrometric reference 
frame, based on specially-designed observations obtained with the two main imaging instruments 
ACS/WFC and WFC3/UVIS onboard the Hubble Space Telescope (HST). Internal comparisons indi¬ 
cate a sub milli-arc-second (mas) precision in the positions and the presence of semi-periodic system- 
atics with a mean amplitude of ~0.8 mas. We combined these observations with numerous archival 
images taken with WFPC2 and spanning 17 years. The precision of the resulting proper motions for 
well-measured stars around the massive cluster R136 can be as good as ~20 /uas yr _1 , although the 
true accuracy of proper motions is generally lower due to the residual systematic errors. The observed 
proper-motion dispersion for our highest-quality measurements is ~0.1 mas yr . Our catalog of 
positions and proper motions contains 86,590 stars down to V ~ 25 and over a total area of ~70 
square arcmin. We examined the proper motions of 105 relatively bright stars and identified a total of 
6 candidate runaway stars. We are able to tentatively confirm the runaway status of star VFTS 285, 
consistent with the findings from line-of-sight velocities, and to show that this star has likely been 
ejected from R136. This study demonstrates that with HST it is now possible to reliably measure 
proper motions of individual stars in the nearest dwarf galaxies such as the LMC. 

Subject headings: astrometry - galaxies: Magellanic Clouds: individual (30 Dor) 


1. INTRODUCTION 

A number of proper-motion measurements are avail¬ 
able for the Small and Large Magellanic Clouds - the 
Milky Way’s nearest irregular dwarf galaxies - which 
address the ir bulk absolute proper motion and ro¬ 
tation (e.g. iKallivavalil et al.1 1 201 .'I lYieira et al.l 120101 
Ivan der Marel fe Kallivavalill 120 1 *lb The accuracy of 
relative proper motions for individual stars with HST 
has been steadily increasing and now is approaching 
the ^0.05 mas yr -1 level, de monstrated fo r a sample 
of Galactic globular clusters (iBellini et al.l 120141 ). This 
opens up an opportunity to measure reliable individ¬ 
ual proper motions for fast-moving stars in the LMC. 
One such class of objects in the Magellanic Clouds are 
massive OB stars, some of which are thought to move 
at velocities close or exceeding 100 km s _1 and, hence, 
should be detectable from HST observations spanning 
two decades. These short-lived massive luminous stars 
play a major role in galaxy evolution by affecting the in¬ 
terstellar medium via chemical enrichment and other pro¬ 
cesses. One of the best sites for detailed studies of such 
stars is the 30 Doradus (hereafter - 30 Dor) region of the 
LMC, which harbors a rich population of O-type stars, 
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both in the form of star clusters and isolated objects 
( Walborn fc Blades 1997 b The existence of the latter is 
puzzling because it contradicts the governing paradigm 
that stars are born in clusters. An additional poorly- 
understood feature of O-type stars is the very high fre¬ 
quency of bi narity and higher- o rder multiplicity among 
these stars (|Mason et al.l 120091: iSana et all 12012th The 
30 Dor region is unique in the sense that it is the near¬ 
est extragalactic site where such stars are concentrated 
and yet are bright enough to permit high-resolution spec¬ 
troscopy with ground-based facilities. Taking advantage 
of this favorable situatio n , the VL T-FLAMES Tarantula 
Survey (VFTS; I Evans et al.l l2011h has obtained multi¬ 
epoch spectroscopy o f over 800 massive OB stars. A 
series of papers (c.f.. iSanalI2013T) based upon these ob¬ 
servations have addressed such questions as: isolated 
high-mass-star formation, internal-velocity dispersion of 
a massive cluster, rotational characteristics of OB stars, 
the fraction of binaries and triples among OB stars. 

One of the most interesting phenomena among OB- 
type stars is the presence of runaway stars (stars with 
large peculiar velocities relative to their general velocity 
distribution). The space frequency of such stars in the 
Milky Way has been studied bv iStond ( 1991 b showing 
that more than ~40% of O-type stars could be runaways, 
while for the B-type stars this fraction drops tenfold. 
Two basic scenarios have been proposed to explain the 
presence of OB runaways. Assuming that massive stars 
are born in star clusters and knowing that mo re tha n 
half of Galactic O-typ e stars are binaries (|Sana et al.1 
l2012tfChini et al1l2012l b it is thought that, if one compo¬ 
nent of a binary explodes as a core-collapse supernova, 
then the remaining star may attain a large kick veloc- 































2 


PLATAIS ET AL. 


ity and, hence, be ejected fr om a cluster (e.g., iBlaauwl 
I196H iHooeerwerf et al.1 l2000f) . The second scenario as¬ 
sumes that a high-mass star is likely to be ejected from a 
cluster vi a the dynamical me chanis m of 3-body encoun¬ 
ters le.g.. lGvaramadze fc Gualandrisil2011 h If a massive 
young star is found in isolation, then there are two pos¬ 
sibilities: it could be a runaway and in such a case it 
should have a parent star cluster not too far away, or it 
could have been formed in isolation and would show no 
kinematic signature of a runaway. This can be tested by 
measuring the line-of-sight (LOS) velocities and proper 
motions. The latter can be translated into velocities, if 
the distance is known with a reasonable accuracy. At 
the distance of the LMC the depth effect is negligible, 
thus making it ideal for such measurements. Further¬ 
more, proper motions have a compelling advantage over 
line-of-sight velocities since they provide the direction of 
motion in addition to its amplitude. In other words, they 
can constrain the probable site of origin of a runaway 
star. 

The 30 Doradus region of the LMC is an ideal envi¬ 
ronment to test models for runaway stars as it has the 
highest density of massive stars in the LMC and includes 
a number of apparently isolated very m assive stars, 
e.g., the confirmed runaway star VFTS 016 (lEvans et al.1 
l20m Indeed VFTS 016 is a prime candidate for the 
dynamical ejection mechanism due to its extreme youth 
and high mass, implying that the possible origin for this 
object is Radcliffe 136 (R136) cluster itself. However, 
this cluster is too young (~2 Myr) to have yet had a 
supernova explosion. Only accurate proper motions can 
constrain the transverse velocity of VFTS 016 and, thus, 
test the assumption that it is a former member of R136. 

2. HST OBSERVATIONS OF 30 DOR 

We started an HST astrometric survey of the 30 Dor 
region (GO-12499; PI: D. Lennon), designed using ex¬ 
tensive prior knowledge about the massive stars in this 
area. This project is part of the HSTPROMO collabora¬ 
tion a set of HST projects aimed at improving our 
understanding of the dynamics of stars, star clusters, 
and galaxies in the nearby Universe through measure¬ 
ment^ of proper motions and their interpretation (e.g., 
Ivan der Marel et all 201 1 i. The HSTPROMO collabora¬ 
tion draws on techniques, tools, and expertise developed 
over two decades of astrometric research with various 
HST imaging instruments. 

The ground-based VFTS provides multi-epoch optical 
spectro scop y of q^800 OB stars in the Large Magellanic 
Cloud (lEvans et all 12011 ). This spectroscopic survey is 
centered around the massive and young star cluster R136 
and it extends outwards up to r ~10'. Our matching 
HST proper-motion survey was optimized to cover one 
quadrant of the VFTS field-of-view (FOV), the entire 
extent of the R136 cluster and its adjacent structures 
of high stellar density, and to include a bona fide run¬ 
away star VFTS 016. Figure [T] shows the spatial cov¬ 
erage of our proper-motion survey with respect to the 
distribution of VFTS targets. The crosshairs indicate 
the location of the most massive Wolf-Rayet star R136- 
al, usually adopted as a center of the entire cluster R136 

5 For details see the HSTPROMO homepage at 
http://www.stsci.edu/~marel/hstpromo.html 


(|Evans et al.ll201ll iHenault-Brunet et al.ll2012f) . 

Our new HST observations of 30 Dor were obtained 
with two imaging instruments: the Wide Field Camera 
3 (WFC3) and its UVIS channel - sensitive at ultraviolet 
and optical wavelengths - and the Advanced Camera for 
Surveys (ACS) through its Wide Field Channel (WFC). 
In order to maximize the field coverage, both instruments 
observed in parallel, with the WFC3/UVIS designated as 
the primary camera. The FOV of WFC3/UVIS covers 
~7 square arcmin of the sky, while the ACS/WFC cov¬ 
ers ~11 square arcmin. The angular separation between 
the centers of these two cameras in the focal-plane is ~6'. 
Such a configuration requires a minimum of two pointings 
to achieve a contiguous, quilt-patterned field coverage by 
both cameras. Formally, the entire VFTS field can be 
covered by ~17 appropriately-spaced pointings, but that 
would result in apparent “holes” in coverage and very 
little overlap between any two adjacent images (frames). 
For the sake of high-precision astrometry with HST, it 
is advantageous to have a ~50% overlap between ad¬ 
joining frames, thus providing a reasonable compromise 
between the minimum number of needed pointings and 
the desire to image the largest possible contiguous area 
of the sky. In the case of the 30 Dor area, this amounted 
to 15 pointings and ~80% of VFTS targets covered, as 
shown by Figs. 11121 Thus, the entire astrometric field in 
its global metachip coordinates is approximately 13' x 16' 
in extent, rotated by ~35° from the East to North. 

In Figure 0 a semi-rectangular pointing pattern 
is generated using the sky location of the primary 
WFC3/UVIS frames. Each of the 15 blue squares rep¬ 
resents a separate pointing. Each adjacent pointing is 
half-a-frame apart in the long direction of a field and it 
has a ~22% overlap in the short direction. Each pointing 
is marked by the visit number (set) in our HST program 
and corresponds to one complete orbit. Sets 04, 05, and 
53 are rotated by ~10° with respect to the bulk of the 
sets. This rotation was necessary in order to have an ap¬ 
propriate pair of guide stars. The parallel observations 
with ACS/WFC (red squares) produced a different pat¬ 
tern of overlaps, where the largest overlap reaches ~43% 
only. 

To enhance the sampling of the point-spread function 
(PSF) and have a better handle on short-scale astromet¬ 
ric distortions, we obtained multiple dithered exposures 
at each pointing, as shown in Fig. [3] We used equal- 
size consecutive steps of 10" along both pixel axes of 
WFC3/UVIS. In this figure, the sub-pointing B defines 
our base pointing for Set 01 shown in Fig. [2j At the 
sub-pointing 1, one short and one long exposure were ob¬ 
tained. The remaining sub-pointings were used to obtain 
a single long exposure only. The applied dither pattern 
for WFC3/UVIS resulted in corresponding shifts almost 
entirely along the y pixel axis for ACS/WFC. The short 
exposures amounted to 35 s for WFC3/UVIS and 32 s 
for ACS/WFC. For WFC3/UVIS, the long exposure was 
either 699 or 507 s, while for WFC3/UVIS, it was either 
640 or 377 s. The observations of 30 Dor were made 
on 2011 October 3-8, with the Set 53 obtained on Oct 
29, 2011. While the timing sequence follows the order of 
observational sets, one set (04) was taken right after set 
16. 

All observations were obtained with the F775W filter 
(similar to the Sloan Digital Sky Survey’s i filter), which 
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was considered the best choice for astrometric purposes 
based on our prior experience. We also chose F775W 
because it provides the best chance of measuring accu¬ 
rate positions for the brightest stars, which are saturated 
even in the short exposures. A number of temporally- 
stable features (“bumps”) near the core of a saturated 
PSF are best defined in this filter, thus providing an op¬ 
portunity to use them for astrometry. Additional details 
about the m apping of the 30 Dor region are described in 
ISabbi et~ahi ( 2013 b That paper also provides an image 
mosaic of the field and preliminary photometric proper¬ 
ties of identified sources. Here, we focus on the initial 
analysis of proper motions based upon comparison to 
the existing HST archival data. As part of GO-12499, 
we have recently obtained a succesful additional epoch 
of astrometric HST observations which will serve to de¬ 
rive proper motions over the entire FOV indicated in 
Fig. [2] These observations and a new set of complete 
proper motions will be discussed in a future paper. 

2.1. Archival HST observations of 30 Dor 

The richest collection of archival HST images for the 
30 Dor region has been provided by the now-retired Wide 
Field Planetary Camera 2 (WFPC2). We selected and 
examined a total of 223 frames obtained in the course of 
13 General Observer (GO) programs during 1994-2003. 
Table |T] provides basic characteristics of these frames, 
such as the GO Program ID, the total number of frames, 
the filter selection, the mean epoch of observations, and 
the target coordinates (in time: hr, min, sec for Right 
Ascension and in angular units: deg, min, sec for Decli¬ 
nation). A program is split up if the target coordinates 
differ significantly, that is more than by a half of the 
PCI’s angular size (the smallest among the four WFPC2 
detectors). For GO 8059 (epoch=2000.25), there are nu¬ 
merous frames with a slightly varying pointing, hence a 
range of coordinates is indicated. There is a large vari¬ 
ety of exposure times, ranging from 3 to 1400 s. This en¬ 
ables us to cover the entire magnitude range of our target 
stars, albeit unevenly across the FOV. The spatial dis¬ 
tribution of the WFPC2 observations is dominated by a 
large number of pointings on the rich star cluster R136. 
The remaining pointings are random, mainly dictated 
by the HST configuration at the time of executed pure- 
parallel observations, such as for GO 8059. We note that 
the largest time baselines with respect to our new data 
reach 17 years, allowing determination of proper motions 
with high precision. 

3. DATA REDUCTION 

In this section, we present a detailed description of 
all the steps required to proceed from fiat-fielded frames 
with three HST imaging instruments to constructing 
astrometric catalogs in the system of a geometrically- 
correct reference frame and, finally, to calculating proper 
motions. At the beginnin g of this pro j ect, w e adopted 
the initial source list from ISabbi et al.l (I2013T) with over 
87,000 stars brighter than the instrumental magnitude^ 

6 Instrumental magnitude is defined as —2.51ogS(DN X Gain), 
where DN is the number of “electrons” counted under the best-fit 
ePSF within the inner 5x5 pixels and the resulting sum is then 
scaled up to represent the total flux within the radius of 10 pixels. 
Note that the ACS/WFC and WFC3/UVIS CCD detectors have 


of —9.3 mag (corresponding to S/N~72), detected in all 
available long-exposure frames. The positions of these 
stars are in the system of WFC3/UVIS pixels, mosaiced 
i nto a master list as described in Sect. 4.1 of ISabbi et al.1 
(12013D . The resulting mosaic is ~26,000 x 20,000 pixels 
large. The first step for us was to start with this master 
list and then construct a reference frame (catalog) with 
higher astrometric accuracy. Our approach provides es¬ 
timates of internal precision at various stages of building 
an astrometric catalog. 

3.1. Image centroids and geometric distortion 

The recent epoch GO-12499 science frames are avail¬ 
able in various output formats from the STScI standard 
data reduction pipeline. We used the bias-subtracted, 
dark-subtracted, flat-fielded, and charge-transfer ineffi¬ 
ciency (CTI) corrected frames, encoded as _flc.fits 
files. All WFC3/UVIS and ACS/WFC frames were con¬ 
verted into a metachip format with all real flux data 
translated into the space-saving inte gers. Th e em pirical 
pixel-based correction for CTI fe.g.. lAnderson &; Bedinl 
I2010D in the WFC3/UVIS data was run as a stand-alone 
option, but this software is expected to be integrated into 
the STScI pipeline soon. 

In order to calculate the centroids and fluxes for 
all detections in ACS/WFC images, we used the 
publicly available softw are code img2xym_WFC.09x10 
^Anderson fe Kingll2006D . modified to apply the neces¬ 
sary corrections for geometric distortion in situ. This 
code employs empirical, spatially-variable point-spread 
function (ePSF) to measure precise position and instru¬ 
mental magnitude for each of the detections (star, ex- 
tragalactic object, or artifact). The ending ”09x10” in¬ 
dicates that there are 9x5 fiducial ePSFs for each of 
the ACS/WFC chips. Because our observations were 
made after the Servicing Mission 4, which restored the 
functionality of ACS/WFC camera, it was decided to 
obtain new ePSFs for the filter F775W based entirely 
on the observations of 30 Dor. A combination of small 
and large dithers among the images was instrumental in 
constructi ng reliable pixel-phase-bia s-free average ePSFs 
detailed in lAnderson fc King! ( 2006 ). We also made im¬ 
provements in the corrections for geometric distortion in 
this filter by updating a look-up table and accounting for 
time-dependent variations in the linear terms (lAndersonl 
120071: lllbeda et al.ll20H . 

For WFC3/UVIS, a sparser array of 7 x 8 fiducial 
ePSFs was adopted because the FOV of this camera 
is smaller and spatial variations of PSF are less pro¬ 
nounced than those for ACS/WFC. The WFC3/UVIS 
library ePSFs were constructed from the observations 
of globular cluster lu Centauri taken shortly after the 
installation of this camera (Anderson, in preparation). 
To correct for geometric distortion in the positions ob¬ 
tained from observations with WFC3/UVIS, we used the 
high -precisio n, t hree-p art correction routine established 
hv dBellini et al.1 (|2011l) . 

Besides the measured pixel positions of all detected 
objects, instrumental magnitudes, and the geometric- 
distortion corrected coordinates, the codes mentioned in 

a much better dynamic range than that of the WFPC2 due to the 
better resolution of analog-to-digital converter: 16 bits versus 12 
bits for WFPC2. 

































4 


PLATAIS ET AL. 


this section deliver an important quality parameter of 
the ePSF fit to each image qfit the sum of abso¬ 
lute differences between the observed pixel values and 
the model, divided by the total flux. It characterizes the 
degree of mis-match between the fitted ePSF and the ac¬ 
tual profile over the central portion of a source. The qf it 
can be as low as 0.01. Fainter sources with qfit>0.5 
appear to be too noisy for accurate astrometry. In addi¬ 
tion, we used a plot of qfit vs. instrumental magnitude 
(see Fig.0 to cull apparent artifacts such as cosmic rays 
by defining an empirical curve which separates genuine 
stars from numerous false detections. Extended sources 
such as external galaxies having intrinsically higher qfit 
values might be lost in this partition, but those are not 
targets of our study. A few surviving contaminants were 
identified in the later cross-idenfications among the over¬ 
lapping frames. 

3.2. Astrometric Reference Frame 

The standard practice in working with partially over¬ 
lapping small images of the sky is to build an over¬ 
arching astrometric reference frame, covering the entire 
FOV. This was achieved by hierarchical accumulation 
of distortion-corrected positions from all long-exposure 
ACS/WFC and WFC3/UVIS frames of GO-12499, start¬ 
ing with subsets of the largest-overlap frames to make ini¬ 
tial tiles and expanding them up to a final link between 
the positions from both HST imaging instruments. In 
all solutions, we used only the well-measured stars with 
instrumental magnitudes in the range —14 < TOF 775 W < 
- 10 . 

First, for each set shown in Fig. [2} we created 15 lo¬ 
cal catalogs (tiles) into the system of each sub-pointing B 
(Fig. [3]) by applying a linear three-term polynomial least- 
squares solution in each coordinate and then averaging 
the common positions. Next, we chose tile 4 and 12 (cor¬ 
responding to the numbered sets in Fig. 0 for ACS/WFC 
and tile 3 and 12 for WFC3/UVIS as starting tiles for 
building a total of four strips, extending from South-West 
to North-East. Then, for ACS/WFC, a strip encom¬ 
passing tile 4 was linked to the system of coordinates of 
tile 12. Similarly, for WFC3/UVIS, a strip encompassing 
tile 12 was translated into the system of coordinates of 
tile 3. Finally, the entire positional catalog of ACS/WFC 
was translated into the system of WFC3/UVIS cata¬ 
log. This sequence of building a composite catalog en¬ 
sures that the WFC3/UVIS tile 3, which covers the cen¬ 
tral part of star cluster R136, is a “touchstone” of our 
global coordinate system. While all listed above least- 
squares solutions were linear, the last transformation into 
the system of WFC3/UVIS catalog required two addi¬ 
tional second-order terms. Without these terms, the 
offset “center-to-edge” between the positions from the 
two instruments reached up to ~0.6 WFC3/UVIS pixel. 
It is suspected that this offset may have its origin in 
the ACS/WFC distortion-corrected coordinates. Lastly, 
the global coordinate system, rooted in tile 3, was ro¬ 
tated around the point corresponding to RA=5 h 38 m 15®0 
and Dec=—69°08'37" so that the rotated pixel axes are 
aligned with the two axes of equatorial coordinates to 
within If5. Our final astrometric reference catalog con¬ 
tains 64,396 objects, each with 1-12 detections depending 
on the number of frame overlaps, with the average of four 
detections per star. 


Table 0 shows the accuracy of various solutions in the 
sense of an rms scatter in the residuals. The first column 
lists the name of the HST imaging instrument. The sec¬ 
ond column indicates the type of area involved in a solu¬ 
tion. We note that a “tile-strip” area for WFC3/UVIS is 
represented by the two adjacent averaged tiles. A strip 
includes at least 7 tiles and is more sensitive to large- 
scale positional errors. Small dithers within a tile make 
a tile solution sensitive to the remaining small-scale ge¬ 
ometric distortion. The mean rms for the X and Y 
coordinates are listed in columns 3-4. The number in 
parentheses gives the rms scatter of these mean values 
in hundredths of a milli-arc-second. The total number 
of least-squares solutions, N so i, is given in column 5 and 
the average number of common stars per solution is pro¬ 
vided in last column. The last line shows the astrometric 
performance of both instruments relative to each other. 
For both instruments, Table 0 indicates that the high¬ 
est accuracy at ^0.6 and ^0.5 mas is achieved in the 
“tile-strip” solutions. This corresponds to ~0.01 pixel. 
Remarkably, solving the ACS/WFC catalog into that of 
WFC3/UVIS produces reasonably high accuracy, which 
is a clear sign that the applied corrections for geomet¬ 
ric distortion in each instrument are as precise as 0.5-1 
mas, with the exception of the noted quadratic differ¬ 
ences. The residuals of the last solution as a function 
of global pixel coordinates are shown in Fig. [5] There 
are small semi-periodic systematics with a typical ampli¬ 
tude of ^0.02 WFC3/UVIS pixel equal to 0.8 mas on the 
sky. Currently, this is state-of-the-art HST differential 
astrometry over scales of ~15 / . Differential astrometry 
over smaller scales (e.g., in the cores of globular clusters) 
requires fewer calibrations and transformations, and can 
therefore achieve even higher accuracies. 

3.2.1. Intricacies of astrometry with WFPC2 

The WFPC2 camera was the first HST imaging in¬ 
strument on which the ePSF concep t was successf ully 
devel oped for undersampled images (|Anderson fe Kina 
[200011 . Contrary to governing wisdom at the time, it was 
shown that undersampled images can provide extremely 
precise stellar positions so long as the constructed PSF 
accurately represents a star’s flux distribution on the de¬ 
tector. It also solved the nagging issue of pixel-phase 
errors, which inevitably appear in the case of an in¬ 
adequate model PSF such as a Gaussian curve (e.g., 
lAnderson & Kingl2000( l. For all WFPC2 _c0f .fits type 
frames, we used the library ePSF in a 3 x 3 c onfig ura¬ 
tion over each of the WFPC2 chips (|Anderson fe Kingl 
2000) implemented in the code dubbed img2meta. This 
software includes accurate correction f or geometri c dis¬ 
tortio n for each of the WFPC2 chips (lAnderson fc King! 
120031 1. The library ePSFs are available for the following 
filters only: F300W, F336W, F439W, F555W, F606W, 
F658N, F675W, and F814W. In the case of a missing 
library ePSF, we used the available nearest-in-central- 
wavelength ePSF. An open issue for the WFPC2 im¬ 
ager is the potential bias in positions due to the time- 
dependent CTI effect , whose photometric calibration is 
provided bv iDolphln (2000 . 2009 ). While the effects of 
CTI in photo metry with WFPC2 ha ve bee n monitored 
for years (e.g.. iGolimowski fc Biretteu l2010h . to date no 
detailed report discussing possible CTI effects on as¬ 
trometry with WFPC2 is available. Given that the CTI 
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mostly affects faint sources and that the VFTS stars of 
primary interest to our study are intrinsicallly bright, we 
chose to neglect these effects for now. An additional ar¬ 
gument in favor of ignoring the CTI in WFPC2 is the lack 
of a perceptable slope in our LMC proper motions as a 
function of magnitude (Sect. EP1) . However, such a slope 
might be mitigated but not eliminated, if the distribu¬ 
tion of HST roll angles for WFPC2 frames is reasonably 
random. 

The software for both ACS/WFC and WFC3/UVIS 
provides the distortion-corrected pixel coordinates of 
sources in metachip format; that is the pixel coordinates 
from each of the two CCD chips are transformed into 
a single system of coordinates, which then can be read¬ 
ily used in building an astrometric reference frame or 
in calculating proper motions. The WFPC2 has four 
separate came ras which are inherently prone to gradual 
displacement (lAnderson fc King! l2003f ). The conclusion 
of these authors is that the interchip separations can¬ 
not be predicted with accuracies better than ~0.2 WF 
pixel. Also, the inner edges of the WFPC2 CCD chips 
are stron gly vignetted and m ust be excluded from all 
datasets (jConzaga et a,l.l[20ToL see Table 1.2). 

We further explored the feasibility of constructing a 
metachip frame by calculating the chip constants (co¬ 
ordinate offsets ex, cy , and rotation angle 0 around 
th e point x=425 y= 425 pix els) using the prescription 
of lPlatais et ah ' (l2002i 2006 ) a nd the astrometric refer¬ 
ence frame, described in Sect. 13.21 Following this algo¬ 
rithm, first we found the best reference stars from least- 
squares solutions of each chip into the reference frame 
using only a linear three-term polynomial for each of the 
axes. Then, these pre-selected reference stars were used 
in a global least-squares solution using only the linear 
terms for each of the axes but this time for all chips 
simultaneously. Prior to this global solution the scales 
of all chips were adjusted to that of chip WF3, which 
is the adopted reference chip for the WFPC2 camera. 
The chip constants were found by iterative sampling in 
the (cx, cy , 0) space for each chip (starting from PCI to 
WF4, WF2, and ending with WF3) and searching for a 
global minimum of x 2 ■ Minimization was stopped when 
the change in cx and cy was less than 0.002 pixel and in 0 
- less than 0"5. To calculate the mean chip constants, we 
used 219 WFPC2 frames ranging in time from 1994.004 
to year 2003.1313. A variable portion of the chip con¬ 
stants A Gapx, A Gapy, ARot with respect to chip WF3 
is shown in Fig. [6] as a function of time. We conclude 
that the angular orientation of the chips is stable over the 
considered period of time. A rather large scatter between 
selected epochs mainly reflects the uncertainties due to 
the varying number of reference stars. However, there are 
apparent changes in relative positions between the chips 
that is a non-linear function of time. There are also some 
smaller ~0.2 pixel fluctuations in the chip separations. 
The pattern of changes and their numerical values are in 
excellent agreement with the results of a similar study by 
lAnderson fc Kind (12003H . It is interesting to note that 
the gap evolution for WF2 and WF4 appears to be identi¬ 
cal (in Fig. [6] compare A Gapx for WF2 with A Gapy for 
WF4), considering the swap of a readout direction. Un¬ 
fortunately, our data are seven years short of the WFPC2 
decommissioning time in 2009. An unpublished study 


(L. Bedin, private communication) indicates that there 
is a turnover in the gap expansion around the year 2004. 
Therefore, the incomplete curves shown in Fig. [6] can¬ 
not be extrapolated beyond the year 2003. In terms of 
proper motions with WFPC2, there are two alternatives: 
1 ) pixel positions from each chip could be considered in¬ 
dependently; or 2) one could construct an astrometric- 
grade global (metachip) coordinate system from all four 
CCD chips using the chip constants, inferred from the 
same frame. In the case of linear transformations that 
are using a three-term polynomial, both approaches pro¬ 
duce similar levels of astrometric accuracy. Considering 
this, we opted for the simpler and more straightforward 
independent-chip approach. 

The multiple archival images of the 30 Dor region ob¬ 
tained with WFPC2 provide an opportunity to examine 
the high-frequency part of the geometric distortion for 
each of the four chips, again with respect to the astromet¬ 
ric reference frame. This is facilitated by the very small 
proper motions of individual stars in the LMC and the 
anticipated tiny positional displacements over a nearly 
two-decade extent of the archival WFPC2 epochs. This 
time we accumulated residuals from linear solutions for 
stars brighter than instrumental magnitude —5.5 and not 
exceeding the averaged residual per each WFPC2 pixel 
by more than 0.1 for PCI and 0.2 pixel for all WF cam¬ 
eras, accordingly. Since the astrometric reference frame 
is on the system of WFC3/UVIS pixels, residuals (and 
the following 2- D maps) are necessarily in the same sys¬ 
tem. Similar to lAnderson fc Kirn] (2003), we averaged 
residuals on a 17 x 17 grid with steps of 47 pixels and 
in bins with a radius of 42 pixels from each gridpoint. If 
a bin contained less than 50 residuals, it was discarded, 
thus resulting in empty bins near the corners of chips 
and in the vignetted parts of chips (Fig. [TJ). The aver¬ 
aged residuals indicate locally-correlated patterns which 
is a hallmark of the residual high-frequency component 
in the geometric distortion. On the system of chip’s na¬ 
tive pixels (as measured and corrected for geometric dis¬ 
tortion), even the largest averaged residuals do not ex¬ 
ceed ~0.03 pix. Therefore, the derived spatial maps of 
residual geometric distortion are useful mainly for studies 
targeting extremely high astrometric precision, such as 
are presented here. We subtracted the inferred residual 
geometric distortion from the WFPC2 pixel coordinates 
prior to calculating proper motions. We note that at this 
level of precision other effects, e.g, the filter-dependent 
part of geometric distortion, might also become signifi¬ 
cant, but we have not explored this. 

Finally, the same residuals from the linear solutions of 
WFPC2 frames into the astrometric reference frame were 
used to obtain the curves of expected standard positional 
error as a function of instrumental magnitude. Due to 
the high similarity of centering errors among all three 
Wide Field cameras of WFPC2, their residuals were com¬ 
bined. The residuals were binned by magnitude and in 
each bin a median was calculated, which then was con¬ 
verted into the estimate of standard deviation by apply¬ 
ing a factor of 1.4825. The resulting curves of expected 
standard error are given in Fig. [8] Although these resid¬ 
uals include proper motions, their contribution appears 
to be insignificant due to the remoteness of LMC and 
the very low expected apparent velocity dispersion in the 
tangential plane. 
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3.3. Proper Motions 

With the benefit of having a comprehensive astromet¬ 
ric reference catalog, generated from our long-exposure 
frames of GO-12499, it is now possible to link any exist¬ 
ing frame in the region of 30 Dor to this catalog. For 
the sake of convenience in finding the common stars, 
we supplemented the cata log with miss ing bright stars 
(m F 775 W < —14) from lSabbi et al.l (120131 1. Our reference 
catalog is in the system of distortion-corrected pixel co¬ 
ordinates of frame ibsf03irq taken with WFC3/UVIS. 
First, we used the global distortion-corrected pixel po¬ 
sitions and rotated them to match closely the RA and 
Dec directions. Second, the J2000 equatorial coordi¬ 
nates were calculated for a total of 67,391 stars down 
to instrumental magnitude of m.F 775 W = —10 that de¬ 
fine our reference catalog. T he correspo nd ing so lutions 
into the UCAC3 positions (iZacharias et al.l[2010f) . using 
a linear 3-term polynomial and a total of 623 common 
reference stars in each coordinate, produced standard de¬ 
viations CTRAcos(Dec) = 0"120 and ODec = 0'.'117. We note 
that in the region of the LMC, the latest UCAC4 cata¬ 
log (jZacharias et al.ll2013fl has no improvement over the 
UCAC3 because the ground-based proper motions of the 
LMC stars are too noisy over small areas such as the 
30 Dor. Then, the fits-header information of each tar¬ 
get frame was used to calculate equatorial coordinates at 
the vertices of the metachip’s FOV. These coordinates al¬ 
lowed us to select a well-defined list of link stars, enabling 
the transformation of a target frame into the system of 
a reference catalog. 

3.3.1. Finding Common Stars 

In order to find the common stars between a chosen 
frame and a list of link stars, we used the so-called vecto¬ 
rial differences in the pixel space drawn from the bright¬ 
est 100-150 stars in both lists: the matching objects pro¬ 
duce parallel vectors that appear as a significant clump 
in the 2-D vector space among the randomly distributed 
false matches. A total of 150 2 =22500 vectors can be 
generated for a sample of the 150 brightest stars, which 
appears to be a good compromise between the required 
computational time and a reliability of the clump’s detec¬ 
tion. The 2-D search for the vectorial clump was substan¬ 
tially sped up by working with sorted arrays of vectors. 
Once a clump of matching stars is identified, the next 
step is to produce an initial least-squares solution (nor¬ 
mally, using a 3-ternr linear polynomial model in each 
coordinate) and then to improve it by including all stars 
with reliable positions. The last set of polynomial coef¬ 
ficients (including all significant higher-order polynomial 
terms) is then used for the final transformation of a tar¬ 
get frame into the system of a reference frame, which is 
a basic step in deriving proper motions. 

The technique of vectorial differences is extremely ro¬ 
bust and works for heavily contaminated samples and/or 
small fractions of overlapping areas. There are, however, 
two preconditions for a success: i) the scale must be the 
same in both systems of pixel coordinates, and ii) the 
axes should be collinear. The first condition is trivial 
to satisfy by applying an appropriate scale factor. The 
second can be less tractable, especially for cases of an 
arbitrary roll angle of HST in combination with a low 
fraction of common stars. For sheer simplicity and con¬ 


clusiveness, we chose to sample all possible angles be¬ 
tween the two coordinate systems, from 0° to 360° at 
0?5 steps (Fig. [9]). This also provides a means to sta¬ 
tistically estimate the significance of a detected clump. 
The software developed for this analysis has two modes: 
a) an interactive mode, enabling a visual check of the 
residual plots and allowing to toggle polynomial terms, 
which is helpful in handling complicated cases; b) an 
unattended automatic mode, to process large numbers 
of frames. The output consists of XY pixel coordinates 
in the system of the astrometric reference catalog and the 
measured apparent magnitude for all objects qualified as 
real detections (see Fig. 0]). 

For this study, there are a total of 149 ACS/WFC and 
WFC3/UVIS first-epoch frames (images), thus yielding 
one catalog for each frame. All combined, these consti¬ 
tute a global master list of 396,616 objects in the area 
of 30 Dor with 14.5 < TOF 775 W < 25.5. Apparent mag¬ 
nitudes in VEGAMAG system are obtained by adding 
a photometric zeropoint listed in HST Data Handbooks 
for each of the instruments. No correction is applied 
for the differences in the aperture size. A star in the 
master list may have a single or up to 17 multiple de¬ 
tections (found within the radius of two WFC3/UVIS 
pixels), which is determined by the frequency of frame 
overlap at the location of this star. Multiple detections 
allow us to calculate the mean position, its standard er¬ 
ror in each coordinate, and to exclude apparent outliers. 
The latter were identified by comparing the estimated 
error at a given instrumental magnitude with the actual 
offset from the mean position. We emphasize that our 
master catalog has a high degree of completeness, with 
the exception of brightest stars which are overexposed on 
all available frames, and some fainter stars located too 
close to a bright star. Such a master catalog is handy 
to quickly identify any star on any image but it is not a 
substitute for an astrometric reference frame. 

A total of 174 selected archival WFPC2 frames were 
similarly translated into the system of the astrometric 
referen ce cata log, to produce 174 catalogs. As indicated 
in Sect. l3.2Tl construction of the latter was done for each 
chip independently, using a linear transformation into 
the reference catalog. The number of available reference 
stars per chip ranges between 6 and ^1000, nearly inde¬ 
pendent of chip selection and despite the fact that the 
PCI covers a factor of four smaller area of the sky than 
the remaining WFPC2 chips. There is only a handful of 
shortest 3-5 s exposures through the broadband WFPC2 
filters, thus limiting the number of available reference 
stars per chip to less than 15 stars. We note that 22% 
of the initially selected WFPC2 frames in Sect. 12. II yield 
unsatisfactory astrometric accuracy and were not used 
in calculation of proper motions. In general, our linear 
transformations appear to be well-constrained, which is 
demonstrated by the fact that the average standard error 
of a solution is small - at the level of 0.09 WFC3/UVIS 
pixel (equivalent to 0.04 WFPC2/WF pixel). 

3.3.2. Calculating Proper Motions 

Once all positions in all exposures are on a system 
of the astrometric reference catalog, selecting epoch ob¬ 
servations for a star and the subsequent calculation of 
its proper motion is straighforward. At each position 
from the global master list, we checked which of the 
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~3.2 million detections are within two WFC3/UVIS pix¬ 
els around this positiorQ. For all identified detections 
around each selected position from the global master list, 
first, a linear unweighted least-squares fit was applied to 
the X and Y positions as a function of time, provided 
the epoch span was larger than 3 years and at least 4 
detections were found. The slope from this fit is the 
proper motion and the standard deviation of the slope 
is the error of proper motion. Although the X-axis in 
the astrometric reference catalog is aligned with Right 
Ascension, it has the opposite direction, thus requiring 
a sign flip for the calculated proper motion in A'. Each 
linear fit was tested for potential outliers and the most 
deviant position was eliminated if its departure from the 
fit exceeded 4<r. The solution was then repeated and the 
residuals again checked for a possible outlier, if the limi¬ 
tation of 4 datapoints had not been reached and/or the 
total number of deviant points is not higher than three. 
In fact, even three rejected datapoints signal that there 
might be an image confusion issue and such a star may 
require a closer examination. 

There was some reasoning behind our decision to ini¬ 
tially use unweighted positions in the fit. First, our 
main targets are intrinsically bright OB spectral-type 
stars which are never under-exposed even on very short 
3-5 s exposures through the F555W and F814W filters 
of WFPC2. Effectively, these stars have the accuracy of 
measured positions dominated by systematics with a neg¬ 
ligible contribution by Poisson noise. This is different for 
faint stars for which shot noise dominates and can be es¬ 
timated. Second, the epoch distribution of observations 
is highly uneven and often nearly all positions at a mid¬ 
dle epoch are slightly displaced by >0.05 WFC/UVIS 
pixel with respect to an overall trend in the positions 
with time. In addition, the spread of datapoints at the 
early epochs sometimes reaches a total range of 0.2-0.3 
WFC3/UVIS pixels. These findings hint to possible sys¬ 
tematics in the WFPC2 frames. Such cases can be iden¬ 
tified by the larger-than-normal proper-motion errors for 
a given number of available observations. 

Despite the robustness of unweighted least-squares fits, 
there are a couple of noteworthy caveats: i) no differ¬ 
entiation is given to the datapoints having clearly dis¬ 
parate positional accuracies, owing to significant differ¬ 
ences in signal-to-noise ratio; ii) the presence of poten¬ 
tial systematics is only loosely indicated by somewhat 
elevated proper motion errors at a given stellar magni¬ 
tude. These issues can be mitigated by “equalizing” the 
datapoints, that is, using their weights. For this, we 
used th e curv es of expected standard errors described 
in Sect. 13.2.11 and given in Fig. [8] The presence of a 
priori error estimates also allows us to calculate y 2 and 
the goodness-of- fit probability Q (subroutine GAMMQ in 
iPress et al.lll992|) . In our calculations of proper motions, 
a low probability Q can be due to underestimated ran¬ 
dom errors of positions and/or the presence of systematic 
errors in the positions. The latter point is illustrated by 
Fig. IT01 which provides the weighted fit to the positions of 
a well-measured brighter star. This star has a total of five 

7 A two-pixel-large displacement over 20 years would correspond 
to — 1000 km s —1 at the distance of the LMC. This restriction may 
end up excluding Galactic stars with relatively high proper motions 
projected onto the 30 Dor region. 


tightly-clumped positional measurements on ACS/WFC 
frames in 2011.7. However, the dispersion of WFPC2 
measurements, obtained with all three WF cameras, is 
higher by nearly an order of magnitude, consistent with 
the notion that at brighter instrumental magnitudes the 
WFPC2/PC1 produces three times less accurate posi¬ 
tions than the ACS/WFC and WFC3/UVIS, but the WF 
cameras - even ~6 times worse. In addition, the obser¬ 
vations with WF2 (epoch=1999.7) in T-coordinate show 
a systematic offset comparing to the observations ob¬ 
tained with WF4 (epoch=2000.7). Despite these system¬ 
atics, the calculated proper motion, px = +0.20 ± 0.02 
p,y = +0.14 ± 0.02 mas yr -1 , appears to be lrigly reli¬ 
able, also indicated by the the goodness-of-fit probability 
exceeding 0.95. 

Conversely, stars with significantly fewer measure¬ 
ments are more susceptible to such systematics and often 
have very low probability Q. We noticed that in the ar¬ 
eas of heavy overlap of ACS/WFC and/or WFC3/UVIS 
frames - as a rule at the edges and corners of chip’s FOV 
(see Fig. [2]) the epoch 2011.7 observations tend to pro¬ 
duce significantly higher dispersion (e.g., Fig. [5]) than the 
predicted standard errors would indicate. This inflates 
the x 2 and, via the probability Q, marks the proper- 
motion calculation as unreliable. In addition, the proper- 
motion errors in such cases can be underestimated. One 
way to mitigate these effects would be to rescale the ex¬ 
pected standard errors of ACS/WFC and WFC3/UVIS 
at brighter magnitudes and optimize the y 2 , especially 
for stars with numerous measurements (e.g., n >25). An 
alternative could be to provide both solutions, with and 
without weights. Normally, with sufficient number of 
observations, such solutions differ insignificantly, mod¬ 
ulo the error estimates. Closing discussion on possible 
systematics, we note that in our WFPC2 frames the so- 
called pixel phase error, due to the variations of ePSF 
over time, is contributing to the error budget at levels 
less than ^0.02 WF PC2/WF pixel, which is in agreement 
with conclusions by lAnderson fc King] ( 2000 1 about this 
effect. 

The final catalog of positions and proper motions is 
based on the weighted fits. The upper limit of acceptable 
proper-motion precision in our catalog of 86,590 stars is 
adopted at a M ~ 1.5 mas yr -1 . As expected, proper mo¬ 
tions derived from fewer detections (n < 9) appear to 
have higher errors. However, nearly a half of the high- 
precision proper motions (<r M < 0.3 mas yr -1 ) are based 
on a smaller number of datapoints (Fig. |TT|). Such cases 
should be considered to be less reliable. We expect to re¬ 
analyze all problematic and low-accuracy-yielding cases 
later with the aid of the GO-13359 second-epoch obser¬ 
vations. 

3.3.3. Exploration of local corrections 

As shown in Fig. [5J residuals from a transforma¬ 
tion of the global ACS/WFC pixel coordinates into the 
WFC3/UVIS global coordinates show small coordinate- 
dependent systematics. Even if the positions from 
WFPC2 observations were perfect, such systematics 
would show up in the proper motions to the extent 
that our linear transformations cannot account for them. 
One common practice to deal with such systemat¬ 
ics is to ap pl y a local ad justm ent to each star (e.g., 
lAnderson fe van der Mare] 120101) . It is expected that 
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there is a characteristic scale-length over which the vari¬ 
ations due to spatial systematics are negligible. For 
ACS/WCS. a 1 00 p ixel box (equal to 5") was adopted by 
lAnderson fc van der Mare] ( 20101 ) as such a scale-length. 
In order to eliminate systematics from proper motions, 
the local stars must be drawn from a population with 
physically consistent bulk motion (e.g., a star cluster), 
so that the differences in proper motions among all used 
local stars are smaller than the sought-after systematics. 
We explored how well this might work for the 30 Dor 
region, where the spatial density of stars varies by or¬ 
ders of magnitude, from the center of R136 out to the 
LMC field stars interspersed with filaments of dust and 
molecular gas. First, we chose only the 11 nearest stars 
with TOF 775 W <23 mag and proper motion error less than 
0.3 mas yr -1 . With these restrictions, the average dis¬ 
tance of the outermost star is ~4", albeit in some cases 
it can be as large as 12". Then, the median proper mo¬ 
tion of these local stars, that is the local correction, was 
subtracted from that of a target star. To evaluate the 
impact of local corrections, we selected a small sample 
(n=50) of VFTS stars with high proper motion acuracy 
(tr M < 0.1 mas yr -1 ) and estimated the dispersion of 
their proper motions before and after the correction. We 
did not find an appreciable decrease in the dispersion 
of locally-corrected proper motions and, therefore, opted 
not to apply these corrections to our proper motions. 

3.4. Poleski et al. Proper Motions 

Recently, iPoleski et al.l (|2012l f published a proper- 
motion catalog of over 6.2 million stars in the direction 
of both Magellanic Clouds. The authors used ground- 
based observations spanning 8 years and obtained ~400- 
700 individual epochs per star. The formal errors of the 
derived proper motions can be as low as 0.2 mas yr -1 , 
thus making them potentially attractive from the stand¬ 
point of potential runaway stars. Nearly 75% of VFTS 
stars are present in this catalog, all with /-magnitudes 
ranging from 13.0 to 17.0 mag. In this magnitude range, 
the Poleski et al proper motions have their estimated 
uncertainties better than 0.3 mas yr -1 . We found 50 
common VFTS stars between Poleski et al. and our 
proper motions with mutually comparable formal uncer¬ 
tainties. If the Poleski et al. error estimates are cor¬ 
rect, we would also expect similar proper-motion dis¬ 
persions, which are dominated by the measurement er¬ 
rors. The distribution of measured line-of-sight velocities 
for the dominant stellar populations in the LMC indi¬ 
cates an upper limit of t he intrinsic velocity dispersion at 
~30 km s -1 (Ivan der Marel et~ahl [200211 . corresponding 
to ~0.1 mas yr -1 at the distance of the LMC. We note 
that the systemic rotation of the LMC has no contribu¬ 
tion to the intrinsic velocity dispersion over the small size 
of the 30 Dor field. A star-by-star comparison (Fig. fl2l) 
shows a 0.16 mas yr -1 dispersion for our proper motions, 
consistent with the estimated proper-motion errors. In 
contrast, the same stars in the Poleski et al catalog show 
a large dispersion at the 0.8-1.3 mas yr -1 level, depend¬ 
ing on whether or not the obvious outliers are eliminated. 
This indicates that their proper motion errors for rela¬ 
tively bright stars are underestimated by a factor of three 
or more. With such hig h act ual mea surem ent errors, the 
proper motions from IPoleski et al.l (|2012f ) appear to be 


inadequate for studies of internal velocities in the Mag¬ 
ellanic Clouds. 

4. DISCUSSION AND APPLICATIONS 

Our main goal is to detect OB runaway stars and iden¬ 
tify a likely site of their origin using our new proper mo¬ 
tions. A more general approach, however, would be to 
identify runaway stars without applying any prior knowl¬ 
edge of photometric, spectral and kinematic character¬ 
istics of stars and then focus on the most interesting 
cases. This approach is not straightforward with the 
current data due to the incomplete coverage of the FOV 
and a significant variation in the accuracy of proper mo¬ 
tions across the field caused by the uneven distribution of 
the first-epoch WFPC2 frames. Therefore, we adopted a 
mixed approach: the statistical properties of our proper 
motions are based on a selection of VFTS stars, which are 
spectroscopically confirmed OB stars, but possible run¬ 
aways are searched among all brighter stars with proper 
motions. 

4.1. Characterizing the global proper motion accuracy 

As indicated in Sect. IQ and Table [U the distribution 
of WFPC2 archival frames is patchy and the combined 
spatial coverage is only ~30% of the size of our astromet¬ 
ric field ('Fig- fUSl) . Our proper motions show no percepti¬ 
ble presence of spatial- or magnitude-related systematic 
errors (Fig. Qj] top panel). The proper-motion distri¬ 
bution narrows significantly if only the higher-precision 
proper motions are plotted fFig. UTl bottom panel). In 
the accuracy range o^ < 0.35 mas yr^ 1 and at m F 775 W < 
18, there are a total of 1262 stars. Among them, a hand¬ 
ful of stars have proper motions exceeding 1 mas yr -1 in 
either coordinate, which we excluded from the analysis. 
At the LMC distance modulus m — M k, 18.50, a proper 
motion of 1 mas yr -1 corresponds to 238 km s . So far 
the studies of LOS velocities for the OB s tars i n dicate 
an upper runaway velocity at ~100 km s -1 (lEvans et alJ 
|2015T> . Hence, it appears unlikely that proper motions 
would reveal runaway stars exceeding 1 mas yr -1 . Most 
likely these fast-moving stars are just the foreground 
Galactic stars. Also, with only two distinct epochs avail¬ 
able for these stars, we are not able to confirm their large 
proper motion yet. 

In our following examination of derived proper mo¬ 
tions we focus entirely on brighter stars (m F775W < 18). 
This sample is scrutinized in two complementary ways: 
i) we estimate the bulk properties of these stars; ii) we 
check the proper-motion fit of selected stars. Finally, we 
searched for possible OB runaway stars and estimated 
the degree of success in finding them. In this magni¬ 
tude range, there is a total of 285 VFTS stars with mea¬ 
sured proper motions. Most of them have formal errors 
less than 0.2 mas yr -1 , however, 38% of these VFTS 
stars have less than 9 datapoints per star. As indicated 
in Sect. 13.3.21 a small number of datapoints may pro¬ 
duce relatively poor proper motions. Therefore, we se¬ 
lected a sample of 109 VFTS stars, requiring that the 
number of datapoints be greater than eight, the proper- 
motion error in both axes is ap < 0.1, and the goodness- 
of-fit probability > 0.1. The vector-point diagram of 
this sample is shown in the top panel of Fig. [15] Then, 
these proper motions, separately in each coordinate, are 
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convolved with a unity Gaussian of the width equal to 
the proper-motion error (never allowing it to be smaller 
than 0.02 mas yr -1 ). These Gaussians are summed up 
and the sum is roughly normalized (bottom panels of 
Fig. USD . The derived smooth distributions are then fit 
with a Gaussian. This provides an important character¬ 
istic of our best proper motions. In both axes, the rms 
scatter around zero proper motion is cr pm =0.10 mas yr -1 . 
This can be converted into the observed velocity disper¬ 
sion. Thus, the observed dispersion of transverse veloc¬ 
ities is ~25 km s -1 . This is a factor of ~3 larger than 
the line-of-sight velocity dispersion derived from single 
VFTS stars (Sana et ah, in preparation), but similar to 
the internal vel ocity d ispersion of the old LMC disk stars 
(Ivan der Marel et all 120021: Ivan der Marel fc Kallivavalill 

I2014T ). While we are not able to resolve the internal ve¬ 
locity dispersion of either 30 Dor or the LMC disk, our 
proper-motion measurement accuracies appear to be suf¬ 
ficient to identify possible runaway stars. 

Finally, a good external check of our proper-motion 
accuracy is provided by the mean measured LOS veloc¬ 
ities of VFTS stars. In order to convert them to the 
equivalent proper motions, first, the global mean LOS 
velocity of 268.5 km s -1 was subtracted from each indi¬ 
vidual velocity. We note that the bulk of OB stars in the 
30 Dor region are mainly concentrated in and around 
R136 (NGC 2070) and the OB association NGC 2060 
(see Fig. |T|). There is a ~10 km s -1 difference in the 
mean LOS between these two major stellar nurseries. 
Therefore, our adopted mean velocity is just an arith¬ 
metic mean and is not representing a particular stellar 
system in 30 Dor. Second, we selected 199 VFTS stars 
with actual proper-motion errors a M < 0.14 mas yr -1 . 
Then, a conversion factor from linear to angular veloc¬ 
ities was applied to the offset LOS velocities from the 
global mean, yielding equivalent proper motions along 
the line-of-sight. For these proper motions, we con¬ 
structed a summed-up Gaussian as in Fig. [15] but this 
time using the corresponding standard error of each star’s 
astrometric proper motion. The resultant Gaussian is 
shown in Fig. [TB] with its width equal to 0.08 mas yr -1 . 
It is smaller only by 20% than the similar widths ob¬ 
tained from our best astrometry shown in Fig. 1151 The 
actual rms dispersion of distribution for the 199 VFTS 
stars from astrometric data alone is 0.12 mas yr -1 . If the 
three-dimensional velocity distribution of the stars were 
roughly isotropic, then one would have expected the two 
distributions shown in Fig. [T6]to be the same. This is not 
the case, although the distributions are reasonably sim¬ 
ilar. The fact that the distribution derived from the ac¬ 
tual proper motions is somewhat broader indicates most 
likely that there are small systematic errors in the proper 
motion data, in addition to random errors. To lowest 
order, the dispersion of the systematic errors can be es¬ 
timated to be < yf (0.12 2 ) — (0.08 2 ) = 0.09 mas yr -1 . 
Hence, there may be systematic proper motion errors in 
this high-quality sample that are comparable to the ran¬ 
dom errors. This needs to be taken into account when 
assessing runaway star candidates. 

4.2. Identifying astrometric high-velocity stars 

If we assume that all stars are born in star clusters 
and associations, then any star with a velocity exceed¬ 


ing the cluster’s escape velocity can be considered as 
a dynamically-ejected “runaway”. However, the short¬ 
lived O stars may escape by another mechanism via 
the binary super nova — which can produce much faster 
runaways. Thus. [Blaauw ( 196 11 for Galactic O stars sug¬ 
gests an empirical lower limit of the excess peculiar ve¬ 
locity at 40 km s -1 . At the distance of LMC this corre¬ 
sponds to a proper motion of ~0.2 mas yr -1 . 

To identify fast-moving stars among the sample of 
brighter stars with tof 775 w < 18, we selected stars with 
a calculated proper motion in either coordinate that dif¬ 
fers from zero at the >5 ct significance level. There are 
105 such stars that also have the goodness-of-fit proba¬ 
bility > 0.2 and the number of datapoints greater than 
eight. The vector-point diagram of these stars is shown 
in Fig. [XT] top panel. A total of 93 of these stars have 
VJ0 photometry from Sabbi et al. (in preparation). A 
VI CMD in the region of 30 Dor is provided in the 
same figure. Note the morphology of the red clump at 
1 < V — I < 2, which is stretched out due to the differ¬ 
ential reddening. 

We can divide the sample in color at V — I = 0.9, 
which approximately separates young blue stars from 
old red stars (70% vs. 30%, respectively). The ob¬ 
served fast-moving red stars may simply be outliers in 
the LMC velocity distribution. The LOS velocity disper - 
sion of old red giants is ~25 km s -1 (|Cole et all [2005'). 
This implies that ~10% of such stars are expected to 
have |Vj > 50 km s -1 in one of the two coordinate di¬ 
rections. Therefore, the observed red stars with high 
proper motions are not of interest in the context of the 
runaway star phenomenon (although they might be of 
intere st for search e s of s tripped SMC stars behind the 
LMC: (olsen et all (j2011l i'l. So, we focus instead on the 
young stars. These have t he intrinsi c ve locity disper¬ 
sion of ~10 km s -1 (Ivan der Marel fe Kallivavalill2014lh 
which makes it highly unlikely to find such stars moving 
in the tail of a Gaussian distribution at \V\ > 50 km s -1 . 

We proceeded with a careful examination of every blue 
and potentially fast-moving star, especially its proper 
motion fit. First, we separated all fast moving stars into 
two groups, each associated with the instrument used for 
the last epoch, that is ACS/WFC or WFC3/UVIS. It 
was suspected that some residual effect of CTI might be 
present in either set of observations. Because of a nearly 
fixed HST orientation during the observations, this would 
readily show up in the vector-point diagram as a distinct 
elongated feature of false proper motions. In fact, there 
is no such feature visible (Fig. [TT] top panel), hence, the 
residual CTI effects can be ruled out. Second, we exam¬ 
ined each proper motion fit and correlated it with a table 
containing detailed information on each datapoint such 
as the filter, exposure time, instrumental magnitude, the 
amount of offset in X and Y, the corresponding chip of 
WFPC2, and the location of the star on this chip. Two 
main culprits for significant systematic errors were iden¬ 
tified: a) the effects of vignetting on WFPC2 chips affect 
astrometry at least up to 100 pixels from the inner chip 
edges. Likewise, the outer edges with pixel coordinates 
>750 might be affected by the geometric modelling ef¬ 
fects (where they reach the maximum) and, possibly, by 

8 The /-magnitude, F775W, in our catalog is ~0.2 mag fainter 
than that calibrated for infinite aperture in Sabbi et al. 
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less accurate corrections of geometric distortion near the 
edge. We, therefore, opted to eliminate all datapoints 
from these ambiguous areas; b) systematic offsets also 
appear to be correlated with a specific chip. We some¬ 
times noticed an unusual scatter if the datapoints are 
drawn from two or more different chips, while their na¬ 
tive pixel coordinates are tighly-clumped on each of these 
chips. Other factors, such as exposure time may have less 
contribution to the systematics, but a use of different fil¬ 
ters usually increases the scatter of datapoints. On one 
occasion, we found three stars with a correlated and en¬ 
larged proper motion in Y, all aligned ~200 pixels away 
from one side of a chip. In addition, if prior to a fit, the 
observations at middle epochs indicate no offset relative 
to the last epoch observations but the earliest WFPC2 
epochs show an obvious offset (hinting at a significant 
proper motion), such a case was marked as unreliable. 
Finally, each candidate star was visually inspected on a 
combined image, constructed from long exposures with 
ACS/WFC and WFC3/UVIS. Only stars with clean and 
well-isolated images were kept for the further considera¬ 
tions. 

After discarding all stars having the noted peculiar¬ 
ities, we were left with a list of only six possible OB 
stars with high proper motions. For most of the high- 
proper-motion stars in Fig. 1171 we have reason to sus¬ 
pect that the ostensibly high proper motion is due to 
some residual systematics. This differs from the case for 
typical stars, for which we found that systematic proper 
motion errors do not generally exceed the random er¬ 
rors (see Fig. USD- Our remaining six candidate runaway 
stars are listed in Table 3. Notes to this table: mag¬ 
nitudes m F 555 W (V) and m f 7 T 5 W (-0 from Sabbi et al. 
(in preparation); Right Ascension (RA) and Declination 
(Dec) in J2000 coordinates; proper motion in each coor¬ 
dinate and its standard error in mas yr -1 ; x' 2 normalized 
to the number of degrees of freedom, that is, Nf ram e — 2; 
the goodness-of-fit probability Q ranging between 0 to 1 
(low Q indicates a non-Gaussian distribution of residuals 
and/or a presence of outliers); and the number of frames 
(epochs) used in the calculation of proper motion. The 
associated weighted least-squares fits of these six stars 
are shown in Fig. [181 

Of the six candidate runaway stars in Table 3, three 
are part of the VFTS survey and the three are not. 
On average, The VFTS stars appear to have somewhat 
larger proper motions, up to 0.4 mas yr -1 in one co¬ 
ordinate, while the non-VFTS stars have smaller ~0.2 
mas yr -1 proper motions. In addition, the VFTS stars 
have only 11-13 datapoints available, while the non- 
VFTS stars have 30-60 datapoints. This kind of pat¬ 
tern might be indicative of some lingering systematics 
for these VFTS stars; possibly, their true proper motions 
could be smaller. In Fig. [151 we marked the candidate 
LOS runaway OB stars among the selected VFTS stars. 
Remarkably, only star VFTS 285 is common between 
the two lists of candidate OB runaways. In light of the 
noted astrometric systematics, we defer a deeper analy¬ 
sis on the runaway star status, statistics, and confidence 
levels to the future study encompassing the most recent 
HST observations. 

Since all six stars are located near the rich star cluster 
R136, we can now explore whether any of them may have 


its origin in this cluster. Figure [T9l shows the location of 
all six OB stars from Table 3. The center of R136 coin¬ 
cides with the zeropoint of the provided gnomonic pro¬ 
jection. For convenience, two dotted circles are drawn 
around this center, indicating the estimates of the half- 
light and the tidal radius. The former is based on a 
single-compo nent K ing mo del fit to t he sur face bright¬ 
ness profile ([McLaughlin fe van der Marell l2005h . while 
the latter is calculated using a two-component fit in 
which only the outer component is a King mqdelj_thus 
providing a better fit at large radius (jMackev & Gilmor^ 
l2003h . All stars but ID 6 = VFTS 285 (d = 2.'36) are lo¬ 
cated within ~1' from R136. One of the most intriguing 
features of the stars in Fig [T9] is the orientation of their 
proper motion vector. Only the likely OB star ID 3 and 
VFTS 285 show proper motions loosely consistent with 
the notion that they might be ejected from R136 with 
VFTS 285 even being outside the estimated tidal radius 
of 217. Star VFTS 285 also has an offset LOS veloc¬ 
ity (—39 ± 4 km s -1 relative to the systemic velocity of 
R136), thus supporting its runaway status. Assuming 
that R136 has zero bulk proper motion (not true but 
estimated to be less than 0.1 mas yr -1 ), our measured 
proper motion of VFTS 285 translates into ~ 110 ± 19 
km s -1 runaway velocity in the tangential plane at a po¬ 
lar angle of <j> = —39° relative to the direction outwards 
the center of R136. In terms of runaway velocity, this 
is consistent with the largest LOS of fset ve lo city ident i- 
fied among the VFTS B-type stars (jEvans et al .1120151) . 
However, due to possible remaining systematics in our 
astrometry, the true proper motion of VFTS 285 might 
be lower. Therefore, we only tentatively conclude that 
VFTS 285 appears to be the first likely O-type (classified 
as 07.5 Vnnn) runaway star in the LMC, detected from 
proper motions. 

In summary, we now have only two reasonable candi¬ 
date OB proper-motion runaway stars, but a more com¬ 
plete picture should emerge after the planned reductions 
of the last-epoch GO-13359 observations with HST. 

5. CONCLUSIONS 

We have studied the 30 Dor region in the LMC with the 
aim of determining individual proper motions of stars. 
This study is a pilot project for deriving relative proper 
motions from HST observations over a significantly larger 
contiguous area of the sky than ever before. Great at¬ 
tention has been given to constructing the backbone of 
our astrometric reductions - the astrometric reference 
frame by joining separate distortion-corrected global 
pixel coordinate systems, obtained from our observa¬ 
tions with two HST imaging instruments, ACS/WFC 
and WFC3/UVIS. We were able to limit internal system¬ 
atics to the level of ~0.02 WFC3/UVIS pixel (0.8 mas) 
over scales of ~15'. This was achieved despite the fact 
that the density of optimal link stars per frame varies 
by a factor of 6 over the 30 Dor field. We note that a 
quilt-like pattern of field coverage with various degrees 
of overlap is a key element in the success of global as¬ 
trometry over a large FOV. 

Next, we explored the potential of using archival 
WFPC2 frames to derive relative proper motions in the 
30 Dor region. A total of 176 WFPC2 frames contributed 
to this effort. Although they cover only one third of our 
total FOV, this provides valuable insights on what kind 
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of proper-motion accuracies can be achieved with a 17- 
year time baseline. Our best proper motions near the 
massive cluster R136 have formal precision as high as 
~20 /xas yr -1 . The true proper-motion accuracy is more 
challenging to quantify, owing to the likely presence of 
residual systematics. Nevertheless, we show that it is 
possible to measure proper motions accurately enough 
to identify fast-moving individual stars in the LMC. 

We derived a proper-motion catalog in the region of 
30 Dor for 86,590 stars. Its most valuable part is the 
higher-accuracy portion, listing ~ 1,300 stars brighter 
than TOF 775 W = 18. This list enabled us to select two 
possible proper-motion runaway stars among a shortlist 
of six OB stars. 

In summary, the main accomplishments of this study 
are: 


1. We provide a detailed description of the construc¬ 
tion of an astrometric reference frame based on one- 
epoclr observations with the HST. Internal compar¬ 
isons indicate a sub-mas accuracy in the positions 
of our astrometric reference frame over scales ex¬ 
ceeding 10'. 


2. We calculated relative proper motions for large 
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numbers of stars down to b ~ 25. This is the 
first higlr-accuracy proper motion catalog in the 
region of 30 Dor, allowing us to search for fast- 
moving stars using proper motions. Our detailed 
analysis of a sample of brighter stars showed that 
the achieved accuracy enables an estimate of indi¬ 
vidual proper motions, if the true relative proper 
motion is larger than ~0.1 mas yr -1 . 

3. There is a hint that one LOS-velocity O-type run¬ 
away star also appears to be a proper-motion run¬ 
away. Star VFTS 285 seems to have its proper 
motion consistent with an origin from the massive 
star cluster R136. This will be verified in the future 
study by adding a final epoch of HST observations 
with WFC3/UVIS and ACS/WFC to the existing 
database. 
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Fig. 2. — Sky coverage obtained with two HST imaging instruments. The Northern semi-rectangular pointing pattern is generated by 
the primary WFC3/UVIS frames (in blue). Each pointing is marked by the visit number in our HST program. The parallel observations 
with ACS/WFC are shown in red. The angular size of entire astrometric field is about IS' x 16 / . 
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Fig. 3.— Dither pattern used for each observational set. Shown is the observational set 01 taken with the WFC3/UVIS camera during 
one HST orbit. The dither step is equal to 10 // along each of the pixel axes. The numbers indicate the sequence of consecutive long 
exposures. The sub-pointing B was adopted as representative for an entire observational set - all other sub-pointings and exposures were 
translated into this frame. 
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Fig. 4. — Separation of star detections from apparent artifacts in a WFC3/UVIS long-exposure frame. All detections above the red 
curve were deleted prior to the following astrometric reductions. The same empirical separation curve was used for both HST imaging 
instruments: WFC3/UVIS and ACS/WFC. 
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Fig. 5.— Positional residuals of stars in the overlap area observed by both, ACS/WFC and WFC3/UVIS. The global ACS/WFC pixel 
coordinates were transformed into the system of WFC3/UVIS global pixels. Small semi-periodic systematics with the amplitude up to 
~0.02 WFC3/UVIS pixels are present. This appears to be the accuracy limit of our differential astrometry over scales of a few times the 
detector size. 
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Fig. 6. — Variations over time in the chip rotation and inter-chip gap size. All chip parameters are relative to the WF3 chip of the 
WFPC2 camera. A zero in the ordinate corresponds to the mean chip parameters derived from a total of 219 WFPC2 frames, taken in the 
area of 30 Doradus. No distinction is made for the filter or exposure time, both of which may add noise t o the datapoints. The observed 
trends or the lack of them are in excellent agreement with the earlier findings by I Anderson fc Kind 1120031) . 
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Fig. 7.— Two-dimensional maps of correlated residuals for WFPC2 chips. These smoothed residuals are on the pixel scale of the 
WFC3/UVIS camera and were applied to the distortion-corrected pixel coordinates of each WFPC2 chip. The size of each vectorial 
residual is magnified by a factor of 500. The maximum size of residuals in each coordinate ranges from ~0.04 pix on PCI to 0.06-0.07 pix 
on wide field (WF) CCD chips. These correct ions can be used to enhance the astrometric precision of the standard distortion solution for 
WFPC2 available in I Anderson fe Kind <120031) . 
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Instrumental magnitude 


Fig. 8. — Expe cted standard cente ring error as a function of instrumental magnitude. The estimates for ACS/WFC and WFC3/UVIS are 
reproduced from IBellini et all 1120141 '). For WFPC2, we derived an error curve for PCI and a single combined error curve for the wide-field 
WF2, WF3, and WF4 cameras. We note that the system of WFPC2 instrumental magnitudes on average is ~2.5 mag fainter than that of 
ACS/WFC and WFC3/UVIS. 
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Fig. 9.— Example of finding a vectorial clump. All possible rotational angles be tween a target frame and the astrometric reference 
catalog were sampled. The position of a peak at 154?5 indicates that angle (see Sect. fT~3.1[) . The full-width-half-maximum of the peak is 
only ~ 1?9. 
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Fig. 10.— Weighted least-squares fit for a 17.9-mag star with a large number of measurements (n=68). Left panel: a fit of residuals in X 
as a function of time; right pane l: the same in Y-coordinate. All residuals are with respect to our astrometric reference frame. A detailed 
discussion is given in Sect. 13.3.21 
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Fig. 11.— Histogram of proper-motion errors. White bins indicate the error distribution for all stars while the dark bins are limited to 
proper motions derived from more than 9 datapoints and cut at cr M < 0.3 mas yr —1 . The distribution of errors is for the proper motions 
in one axis. The vast majority of stars with fewer datapoints are fainter than mp 775 w = 21 and the formal high precision of their proper 
motions seen for some stars might be suspect due to the small number statistics. 
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Fig. 12.— Vector-point diagram of selected VFTS stars. Left panel: a total of 50 VFTS stars in our proper motion catalog with a 
proper-m otion error less that 0.13 mas yr —1 . The proper-motion dispersion in either axis is 0.16 mas yr -1 . Right panel: the same stars 
but from [Poleski et ahl (|2012D . The dispersion of proper motions is at least 0.8-1.3 mas yr —1 , depending upon the trimming of a few stars 
outside the shown limits. Our proper motions are clearly more precise and accurate. 
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Fig. 13.— Distribution of stars with measured proper motions over the FOV. Only every 5-th star is plotted (black dots). The coverage 
pattern is defined by the availability of archival WFPC2 images. The rugged line has the same meaning as in Fig. [T] Large cyan dots 
indicate the location of 105 ostensibly faster-moving stars with proper motions exceeding 5<r in at least one coordinate, where a is the 
actual proper-motion error of a star. Our best sample of possible OB runaway candidates is marked by large red dots. 
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Fig. 14. — Relative proper motions in X of stars in the 30 Dor region as a function of magnitude. The subscripts X and Y are a shorthand 
for RAcos(Dec) and Dec, and are used throughout the text, figures and tables. Upper panel: all 86,590 catalog stars are plotted. Lower 
panel: only stars with cr^ < 0.35 mas yr —1 in both axes and with nf rame > 9 in the proper-motion fit are plotted. There are 37,026 such 
stars. The number of faint stars with unrealistic proper motions (>1 mas yr —1 ) drops dramatically once the low-accuracy fits are deleted. 
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Fig. 15. — Distribution of well-measured proper motions for a total of 109 VFTS stars. Upper panel: a vector-point diagram. The LOS 
candidate runaway OB stars are highligh ted b y the red symbols. Star VFTS 285 is labeled because of its likely status of the LOS and a 
proper-motion runaway O star (see Sect. im Bottom two panels: a sum of unity Gaussians (black points) derived using the measured 
proper motions and their errors. Red curves represent a Gaussian fit to these distributions. In both coordinates the Gaussian width is only 
0.10 mas yr —1 . 
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Fig. 16. — Distribution of proper motions for 199 VFTS stars. The black points show a superposition of unity Gaussians calculated 
from the measured proper motion in X. Similarly, the red curve is constructed by converting each offset-from-the-mean VFTS line-of-sight 
velocity into the equivalent of proper motion, convolved with the independently measured astrometric standard error. The distribution 
derived from the actual proper motions is somewhat broader than that derived from the LOS velocities. This indicates most likely that 
there are small systematic errors in the proper motion data that are comparable to the random errors. 
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Fig. 17. — Possible fast-moving stars. Upper panel: a VPD of 105 bright stars (mp 775 W < 18) with statistically-significant proper 
motions at a 5cr confidence. Red circles show our best six candidate runaway OB stars near R 136; bottom panel: a color-magnitude 
diagram indicating the possible fast-moving stars (large cyan points) and the candidate runaways (red circles). The black points represent 
a selection of stars with I < 20 in the 30 Dor region. The vertical line at V — I = 0.9 approximately separates the domain of OB stars 
from the brighter field RGB stars. 
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Fig. 18, 


Weighted least-squares proper-motion fit for six candidate runaway OB stars. See Sect. 14721 Fig. 1101 and Table 3 for details. 
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Fig. 19. — Location and proper-motion vector of six candidate runaway OB stars. The center of star cluster R136 is at the zeropoint of 
gnomonic projection. To aid the eye, a line showing a possible path of ejection from R136 is drawn for each star. Proper-motion vectors 
are shown by the green arrows. In this plot, V corresponds to a proper motion of 1.5 mas yr —1 . Out of six stars, only two have a direction 
of proper motion that suggests ejection from R136. Dotted circles indicate the estimates of the half-light and the tidal radius. 
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TABLE 1 

WFPC2 OBSERVATIONS OF 30 DOR 


Program ID 

Nobs 

Filters 

Epoch 

RA (J2000) 

Dec (J2000) 

5114 

48 

F555W,F814W 

1994.74 

5:38:43 

-69:06:05 

5584 

4 

F656N,F675W 

1995.50 

5:37:01 

-69:07:14 

5589 

21 

F336W,F555W,F814W,F547M,F502N,F656N,F673N 

1994.01 

5:38:42 

-69:06:00 

6122 

13 

F555W,F814W,F502N,F656N,F673N 

1996.96 

5:38:17 

-69:04:02 

6251 

1 

F814W 

1995.50 

5:37:01 

-69:07:14 

6540 

8 

F656N 

1998.24 

5:38:52 

-69:08:07 

7786 

1 

F606W 

1997.79 

5:37:45 

-69:10:59 

8059 

4 

F606W,F814W 

1999.79 

5:39:08 

-69:10:40 

8059 

4 

F606W,F814W 

1999.80 

5:38:52 

-69:11:08 

8059 

44 a 

F606W,F814W 

2000.25 

5:38:16 

-69:01:24 

8059 

23 

F606W,F814W,F656N 

2001.02 

5:37:48 

-69:08:20 

8090 

8 

F606W 

2000.25 

5:38:11 

-69:01:34 

8163 

8 

F502N,F656N,F673N 

2000.73 

5:38:46 

-69:05:10 

8163 

8 

F502N,F656N,F673N 

2000.31 

5:38:48 

-69:04:40 

8163 

8 

F502N,F656N,F673N 

2000.24 

5:38:34 

-69:06:02 

8163 

6 

F555W,F814W 

1999.76 

5:38:44 

-69:07:44 

8163 

6 

F555W,F814W 

1999.76 

5:38:33 

-69:04:33 

8436 

1 

F606W 

2000.22 

5:38:02 

-69:02:14 

8883 

1 

F606W 

2001.02 

5:37:49 

-69:08:18 

9676 

6 

F606W 

2003.00 

5:38:44 

-69:08:27 


a Mean coordinates given. The actual range in RA is from 5:38:08 to 5:38:24 and in Dec is from —69:01:03 to —69:01:46. 
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TABLE 2 

RMS SCATTER OF VARIOUS REFERENCE-FRAME SOLUTIONS. 


Camera 

Field type 

rms in X (mas) 

rms in Y (mas) 

N so i 

N stars 

ACS/WFC 

frame-tile 

0.66(4) 

0.96(4) 

44 

3200 

ACS/WFC 

tile-strip 

0.58(8) 

0.65(7) 

13 

2800 

ACS/WFC 

strip-strip 

0.83 

0.83 

1 

7903 

WFC3/UVIS 

frame-tile 

0.60(4) 

0.74(5) 

44 

1100 

WFC3/UVIS 

tile-strip 

0.45(2) 

0.52(5) 

13 

940 

WFC3/UVIS 

strip-strip 

0.60 

0.88 

1 

2596 

ACS-WFC3 

strip-strip 

0.56 

0.68 

1 

3322 



Table 3. Candidate OB proper-motion runaway stars (see Section 14.21 for units and details) 


ID 

VFTS 

ro-F775W 

7ttF555W — m F775W 

RA (deg) 

Dec (deg) 

Mx 

MY 

°Vx 


Xx 

X 2 r 

Qx 

Qy 

Nf r ame 

1 


16.836 

0.255 

84.6963497 

-69.1091730 

0.03 

0.19 

0.02 

0.03 

0.32 

0.91 

1.00 

0.67 

62 

2 

350 

14.718 

0.316 

84.6169242 

-69.1080834 

0.37 

0.36 

0.05 

0.06 

1.03 

1.07 

0.41 

0.38 

11 

3 


15.508 

0.308 

84.6776254 

-69.1025824 

-0.02 

-0.18 

0.02 

0.02 

0.32 

0.22 

1.00 

1.00 

42 

4 


17.024 

0.264 

84.7270015 

-69.0998315 

-0.06 

0.22 

0.02 

0.02 

0.23 

0.23 

1.00 

1.00 

31 

5 

373 

15.096 

0.225 

84.6301239 

-69.0980782 

0.26 

-0.21 

0.06 

0.04 

1.29 

0.49 

0.22 

0.91 

13 

6 

285 

15.482 

0.153 

84.5722167 

-69.0950254 

-0.40 

-0.23 

0.06 

0.05 

0.79 

0.56 

0.64 

0.85 

12 


CO 

CO 


Proper motions in 30 Doradus 






